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Probabilistic quantum state transformations can be characterized by the degree of state separation 
they provide. This, in turn, sets limits on the success rate of these transformations. We consider 
optimum state separation of two known pure states in the general case where the known states have 
arbitrary a priori probabilities. The problem is formulated from a geometric perspective and shown 
to be equivalent to the problem of finding tangent curves within two families of conics that represent 
the unitarity constraints and the objective functions to be optimized, respectively. We present the 
corresponding analytical solutions in various forms. In the limit of perfect state separation, which is 
equivalent to unambiguous state discrimination, the solution exhibits a phenomenon analogous to a 
second order symmetry breaking phase transition. We also propose a linear optics implementation 
of separation which is based on the dual rail representation of qubits and single-photon multiport 
interferometry. 
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I. INTRODUCTION 

Quantum information processing deals with changes in 
the state of a quantum system and what they amount to 
in terms of the information encoded in the initial state 
and the transformed or final state. Not any transforma¬ 
tion from a given state, or set of states, to another is 
allowed by quantum mechanics, which sets strong limita¬ 
tions to the processing of information in quantum compu¬ 
tation and quantum communication [T]. Even so, quan¬ 
tum information processing is expected to outperform 
its classical counterpart [2j. Some expectations are al¬ 
ready materializing in quantum cryptography 011] and 
quantum simulation 0 0: and much more is to come 
as experimentalists make progress overcoming decoher¬ 
ence and other issues involved in the implementation of 
quantum technologies. 

Since the evolution of quantum states is central in 
quantum information, we are urged to investigate the 
ultimate limits imposed by nature on state transforma¬ 
tions. In this sense, it has been recognized that prob¬ 
abilistic processing can offer significant advantages over 
deterministic processing. The simplest example is ar¬ 
guably unambiguous discrimination [ZHna, which en¬ 
ables error-free identification of non-orthogonal quantum 
states, provided they are linearly independent m (El- 
More recently, perfect cloning has been proved possible 
when prior knowledge about the possible preparations of 
the state to be cloned is given Cana. In all these cases, 
the price to pay is, of course, that the processing fails 
some times. However, we have means to know that the 
process has failed and we can compute the failure rate. 

The examples above, as well as more recent devel¬ 
opments in quantum replication | T6] and probabilistic 
metrology m, may be just the tip of the iceberg pointing 


at new directions in quantum processing. In this paper 
we will focus on the simplest case of transformations over 
pure states belonging to a given two-state family. Such 
families are characterized by the overlap of their states 
and their prior probabilities. Any transformation acting 
on them gives a new two-state family also characterized 
by the overlap of the transformed states. Whether or 
not the transformation is possible with some given fail¬ 
ure rate depends solely on the value of these overlaps 
and the prior probabilities of each of the states of the 
original family. For zero failure rate, i.e., for determin¬ 
istic processes the final overlap is necessarily larger or 
equal to the initial overlap. However, if some non-zero 
failure rate is allowed, the final overlap can be smaller 
than the initial overlap, and we say that the transforma¬ 
tion increases the degree of separation [ISJ of the original 
states, since the final states become more easily distin¬ 
guishable. Two states become fully separated under a 
transformation if the corresponding transformed states 
are orthogonal, i.e., they have zero overlap. Full sep¬ 
aration is equivalent to unambiguous discrimination in 
the following sense. If the transformed states are orthog¬ 
onal, they can be discriminated with no ambiguity by 
a projective measurement along the rays defined by the 
transformed states. So, the transformation can be used 
to implement unambiguous discrimination with the very 
same failure rate. Conversely, if two states can be un¬ 
ambiguously discriminated, upon identification we can 
prepare any state, in particular a state out of a pair with 
zero overlap. This shows that unambiguous discrimina¬ 
tion followed by state preparation implement any trans¬ 
formation that fully separates two states. So, there is a 
measure-and-prepare protocol that implements any such 
transformation. 

In intermediate situations, where some degree of sep- 
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aration is attained, there are several questions that we 
should answer. If the degree of separation is given, what 
is the optimal protocol, i.e., the protocol that has the 
smallest failure rate? If the failure rate cannot exceed a 
given value, what is the maximum degree of separation 
a transformation can attain as a function of the origi¬ 
nal overlap? And finally, what is the tradeoff between 
degree of separation and failure rate for a given initial 
overlap? These three questions are not independent, of 
course, but because of the impossibility to find a fully 
explicit solution to separation, which would involve solv¬ 
ing sixth degree polynomial equations, each one of them 
must be addressed separately. We provide the answers, 
i.e., the plots of the quantities relevant to each situation, 
in a simple parametric form. This gives a full account 
of the separation problem. The geometric approach, de¬ 
veloped in ib-ii and [ 21 ], proves equally powerful here. It 
encompasses the entire physics in a simple intuitive pic¬ 
ture and lends itself to analytical or numerical studies for 
which it provides a visual guidance. 

A phenomenon analogous to a second order symmetry 
braking phase transition arises in the limit of full sepa¬ 
ration, i.e., when the overlap of the transformed states 
vanishes. This was already noticed in our recent let¬ 
ter m on perfect cloning, which is a particular instance 
of separation since the overlap of the perfect clones is 
necessarily smaller than the overlap of the states to be 
cloned. There, we showed that the failure probability as 
a function of the prior probabilities is an analytic func¬ 
tion if a finite number of clones are produced, but its 
second derivative becomes discontinuous in the limit of 
infinitely many clones. In this limit full separation takes 
place, since the overlap of the clones approaches zero ex¬ 
ponentially as the number of clones increases. Here we 
show that such phase transition is a general feature of 
separation. 

The paper is organized as follows. In Sec. [H] we intro¬ 
duce the separation problem and our notation. We also 


show in detail that a unique solution exists. In Sec. Ill 


we derive the minimum failure probability for a fixed de¬ 
gree of separation as a function of the prior probabilities. 
In the particular case of perfect cloning we recover the 
results of our previous work in m- In Sec. IWl we de¬ 
rive the maximum separation for a fixed failure rate as a 
function of the initial overlap. In Sec. [V] we obtain the 
tradeoff curve between degree of separation and allowed 
failure probability. In Sec. |VI| we provide a physical im¬ 
plementation based on single-photon, multiport interfer¬ 
ometry employing the rail representation of qubits. We 
close with a brief discussion of our results in Sec. rvTTl 


II. SETUP FOR QUANTUM STATE 
SEPARATION 

We can always imagine that a probabilistic quantum 
transformation is carried out by a machine with an input 
port, an output port and two flags that herald the success 


or failure of the transformation. The input \ipi), i = 1 , 2 
is fed through the input port for processing. In case of 
success, states |t/>'), with the desired degree of separation, 
are delivered through the output port with conditioned 
probability pi. Otherwise, the output is in a failure state. 
Conditioned on the input state being \ipi), the failure 
probability is qi = 1 — p^. 

We address optimality from a Bayesian viewpoint that 
assumes the states to be transformed are given with some 
a priori probabilities r)i and 772, 771+772 = 1 - Then a 
natural cost function for our probabilistic machines is 
given by the average failure probability 

Q = 7?i <71 +772172- ( 1 ) 

If | t/ji) and the corresponding transformed states |t/') are 
given, the optimal machine is one that minimizes the 
cost function Q. In this case our aim is to find that 
optimal machine and the minimum average failure prob¬ 
ability Qmin for arbitrary priors 771 and 772 ■ 

A different way of approaching optimality may con¬ 
sist in finding the machine (or machines) that achieves 
the highest degree of separation, namely, minimizes the 
overlap s' := |( 7 /A|t/ 4 )I for given initial states | ipi), sub¬ 
ject to the condition that the average probability Q 
does not exceed some given value, Q max . In this case 
we could further assume that either the initial overlap 
s := Kt/AIt/^)! is given, in which case one can compute 
the tradeoff curve s( n i n (Qmax), or else assume that Q max 
is fixed and compute the curve s+ n (s). It is easy to see 
that s' min (Q max ) and Q max (s' min ) are just inverses of each 
other. 

Whether we approach optimality one way or another 
depends merely on the problem at hand. Hence, e.g., for 
perfect cloning from one initial copy of either | tp\) or 1 7 /> 2 ) 
to n final copies (i.e., |^') = [t/A)®”), the former approach 
is most suitable since the final overlap is fixed, s' = s n , 
and so is the degree of separation attained by the cloner. 
So, in [T 5 ] the solution was given in terms of Q m i n as a 
function of the prior probability 771. However, one may 
need to know what is the maximum number of clones that 
can be produced if the failure rate cannot exceed Q max , 
in which case one takes the latter approach, and com¬ 
putes 77 max = log [s' (Qrnax)] /logs. 

The machine that carries the probabilistic transforma¬ 
tion is usually described by two Kraus operators A succ , 
Af a ii, so that A / ucc d succ + d ) a i 1 A 1 = U [1 , 18 . We 
can think of d succ and df a n as measurement operators. 
The transformation is successfully applied if the out¬ 
come of such (generalized) measurement is “succ”, and 
fails otherwise. Neumark’s theorem provides an alter¬ 
native approach that turns out to be more convenient 
for our analysis. Additional details on this method can 
be found in m- In this formulation, the Hilbert space 
df of the original states is supplemented with an ancil¬ 
lary space extra <S> Mf that accommodates both the 
required extra-dimensions (if necessary) as well as the 
success/failure flags. Then, a unitary transformation U 
(time evolution) from K g df extra g TCp onto !K' g 3 ~Cp is 
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defined through [T 31 HU [H] 


U 1-01)10) = V / PilV ,, i)|ai) + v / 9il^)l a o) ) (2) 

U\if>2)\0) = V^ 2 >|a 2 > + n/9 2 I^)I“o). (3) 

Here the ancillas are initialized in a reference state |0). 
The states of the flag associated with successful trans¬ 
formation |a,) are constrained to be orthogonal to the 
state |ao) that signals failure. Upon performing a pro¬ 
jective measurement on the flag space 3~Cp, the final state 
delivered through the output port of our probabilistic 
machine is either |i/ 4 ), in case of success, or |</>) in case 
of failure. So, the outcome of this measurement tells us 
if the machine has succeeded or failed in delivering the 
right transformed state. On general grounds, optimal¬ 
ity requires |ai) = |a 2 )- Here we choose to consider 
a more general setup where these two states are differ¬ 
ent to include state discrimination, for which the success 
flag states must be fully distinguishable, so (ai|a 2 ) = 0 . 
Likewise, we could consider an even more general setup 
with two failure states \(j)\) and |</> 2 ) in Eqs. © and ([3]). 
This is necessarily sub-optimal since we could proba¬ 
bilistically determine whether we received \if\) or |^ 2 ) 
by applying unambiguous discrimination to the failure 
states | (f>i). Sometimes we would be certain of the in¬ 
put state, in which case we could prepare or \i/j 2 ) 
accordingly, thereby increasing the overall success rate. 

Taking the inner product of Eqs. © and ([3]) with 
themselves shows that our probabilities are normalized: 
Pi T Qi = 1. Similarly, by taking the product of Eq. ([ 2 ]) 
with Eq. ©, we find the unitarity constraint, 

s = \JP\Vi fi + y/qiq ?, (4) 

where (3 = s'|(ai|a 2 )|. Without any loss of general¬ 
ity, in deriving Eq. Q we have chosen (ipilfa), (V4 Wi) 
and (ai|a 2 ) to be real and positive. We note that 
0 < (3 < s, and (3 = 0 for both full separation (s' = 0) and 
unambiguous discrimination (|(ai|a 2 )| = 0 ), whereas for 
optimal separation |(ai|a 2 }| = 1- If Eq. Q is satisfied, it 
is not hard to prove that U has a unitary extension on the 
whole Hilbert space and the Kraus operators, A succ , -Afaii, 
can be obtained by tracing out the ancillary degrees of 
freedom. 

Geometrically, Eq. 0 defines a straight line in the 
< 71 -c /2 plane for fix values of Q and the priors. Using 
Pi = 1 - qi, Eq. © defines curves in the same plane 
characterized by the values of s and (3. In Fig. [T] we 
display these lines and curves for representative values 
of the parameters. For convenient referencing, we gather 
in a lemma all the features of these curves that we will 
need. Points (a)-(d) are straightforward, so only (e) and 
(f) are proven below. 



FIG. 1. Unitarity curves in Eq. pll and the associated sets S /3 
in Eq. |5} for /3 = 0.45 (solid/light gray), j3 = 0.30 (dot¬ 
ted/medium gray), and /3 = 0 (dashed/dark gray). The fig¬ 
ure also shows the optimal straight segment Q = rjiqi + 77292 
and its normal vector ( 771 , 772 ). Plotted for s = 0.6, 771 = 0.17, 
772 = 0.83 and Q = 0.24. 

Lemma 1 (a) For fixed s, Eq. 0 defines a class of 
smooth curves on the unit square 0 < qt < 1 (e.g., solid, 
dashed or dotted curves in Fig. Qp. (b) All these curves 
meet at their endpoints, (l,s 2 ) and (s 2 ,l). (c) At the 

endpoints the curves become tangent to the vertical and 
horizontal lines qi = 1 and <72 = 1 respectively, provided f3 
is not zero, (d) For (3 = 0 the curve is an arc of the hy¬ 
perbola q±q 2 = s 2 (dashed line in Fig. ©). (e) Each of 
these curves and the segments joining their end points 
with the vertex (1,1) enclose the sets (any of the gray 
regions in Fig. [7]) 

<S/3 = {(< 7 i,g 2 )e[ 0 ,l]x[ 0 ,l]: \JP\P 2 P + V 9 K 22 - s > 0 }. (5) 

They satisfy Sp C Sp> if ft < fi'. (f) Moreover, the 

sets Sf 3 are convex. 

Proof, (e) The curve Q is readily seen to be part of the 
boundary of Sp. Assume that /3 < (fi and (gi,q 2 ) £ Sp. 
Then 

y/piP2 P' + V9W2 - s> sjp\pi (3 + ^9192 - s > 0, (6) 

and thus (<7i,g 2 ) £ Sp' ■ (f) To prove convexity let us 
assume that (< 71 ,( 72 ) and (gj,^) belong to Sp. We define 
qi = Mi + (1 — A )q(, where 0 < A < 1. It follows that 
pi := 1 - % = A pi + (1 - A )p(. Since f(x, y) = v fxy is a 
concave function in the unit square {(a;, y) \ 0 < x, y < 1}, 
we have \fq\q 2 > A A /gig 2 + (l — an d> since /3 > 0, 

VP 1 P 2 P > Asjpxp 2 + (1 - A ) s /p' 1 p' 2 fi. Then 


VP1P2 fi + Vqi02 - s > A (^pip 2 fi + sjqfiqi - s) + (1 - A) (VpiP 2/3 + VdW2 - s ) > 0 . 


(7) 
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Thus, ( 91,92 ) € Sp, which proves the convexity of Sp 
for /3 > 0. ■ 

Now that we have characterized the geometry of the 
unitarity constraint, a geometrical picture of the opti¬ 
mization problem emerges (See Fig. [TJ. Eq. |l]) defines a 
straight segment on the square 0 < qi < 1 with a normal 
vector in the first quadrant parallel to (771 , 772 )- For fixed 
a priori probabilities, the average failure probability Q 
is proportional to the distance from this segment to the 
origin ( 0 , 0 ). The intersection of such a straight segment 
with the boundary of Sp provides an admissible unitary 
transformation U and its corresponding failure probabil¬ 
ity Q. Since Sp is convex and the stretch of its bound¬ 
ary given by Eq. Q is smooth, the optimal transforma¬ 
tion, for which Q is minimal, is defined by the unique 
point ( 91 , 92 ) °f tangency with the segment 0 that ex¬ 
ists for any value of the priors and for /? > 0 . So, this 
tangency point determines the minimum failure proba¬ 
bility Q m i n and defines the optimal separation strategy 
through Eqs. p]) and §. 

For /3 = 0 (full separation/unambiguous discrimina¬ 
tion), the right hand side of Eq. Q describes a hyper¬ 
bola for a fixed value of s, 92 = s 2 /qi, corresponding 
to a dashed line in Fig. [lj Its slope, q 2 = —s 2 /9 2 , 
is in the range [— s -2 ,— s 2 ]. A unique point of tan¬ 
gency with the line (|T|) can only exists if the slope of 
this line, — 771 /772, is within this same range, namely if 
s 2 /(l + s 2 ) < 771 < 1/(1 + s 2 ). The tangency point 
is then seen to be (91,92) = ^mm s (77 1 ,77^" 1 ). This 
leads to a minimum average failure probability given by 
Qmin = Qud := 2 \/mV is, where the subscript UD stands 
for unambiguous discrimination. If the slope is outside 
the range tangency is not possible, and then the optimal 
line merely touches the end points of the hyperbola. For 
Pi < s 2 /( 1 + s 2 ), the straight segment 0 pivots on the 
lower end point, (l,s 2 ), as we vary 771 and we have the 
minimum average failure probability as Qud = rji+^s 2 . 
Likewise, for 771 > 1/(1 + s 2 ), the pivoting point is the 
upper end point of the hyperbola, (s 2 ,l), which leads 
to Qud = 771 s 2 + 772- The above can be summarized by 


ty/Vim s , 


1 + s 2 


<Vi< 


1 + s 2 


Q UD — < 


m + S 2 772, 0 < 771 < 


1 + s 2 ’ 


( 8 ) 


ms + m, 


1 


1 + s* 


<m< 1 . 


This expression reproduces the optimal average failure 
probability for unambiguous discrimination [TO] , as it 
should. 

Furthermore, we note that for the second (third) line 
in ([8]) we have pi = 1 — 91 = 0 (j?2 = 1 — 92 = 0), which 
leads to a 2-outcome projective measurement, as only 
the success flag state l^) (|ai)) is needed in Eqs. 0 
and ([ 3 ]). The solution in the first line of Eq. Q is mani¬ 
festly symmetric under the exchange of the input states, 
i.e., under 771 0 772. However, this symmetry is lost in 
the other lines. Instead, the effect of swapping the states 


turns the solution in the second line of Eq. Q into the 
solution in the third line. One can also check that Qud 
is a twice differentiable function of 771 (or 772), with a 
second derivative discontinuous at m = s 2 /(l + s 2 ) and 
771 = 1/(1 + s 2 ). Our geometrical approach shows that 
the average failure probability Q m i n is an infinitely dif¬ 
ferentiable function of 771 for (3 > 0 , since according to 
our lemma, the boundary curve 0 merges smoothly into 
the lines 91 = 1 and 92 = 1. So, it turns out that at 
(3 = 0 a phenomenon similar to a second order symme¬ 
try breaking phase transition takes place. A similar phe¬ 
nomenon was observed in unambiguous discrimination of 
more than two pure states m- 

Our lemma can likewise be used to address optimality 
for given priors m and m and average failure probability 
not exceeding Q m ax , with 0 < Q max < Q ud- First, since 
the unitarity curve is a function of f3 = s' (a^az), we set 
|oq) = |a 2 ), be., (oq|cr 2 ) = 1, to ensure the minimum 
value of s' > 0 for a given j3. Then, it follows from the 
lemma that the minimum final overlap s' > 0 (the maxi¬ 
mum degree of separation attainable), which we call s y mirl , 
is that for which the segment ([I]), with Q = Q max , and 
the boundary of Sp- s > become tangent. Setting the mar¬ 
gin Q max in the range [Qud,1] leads, obviously, to the 
trivial solution s( nin = 0, for such margin would allow 
full separation using unambiguous discrimination with a 
failure rate of exactly Q = Q ud, below the given margin. 

In summary, our lemma provides the solution to op¬ 
timal state separation from a geometrical viewpoint by 
showing that it is a convex optimization problem, for 
which a unique solution exists. Unfortunately, a closed 
form for this solution does not exist for arbitrary prior 
probabilities, since finding the tangency point of the seg¬ 
ment in Eq. ([!]) with the curve in Eq. Q requires solving 
a six degree polynomial equation, as one can easily check. 
In the next sections, we give an analytic solution to state 
separation in parametric form. This solution contains all 
the information one may need in a simple and straightfor¬ 
ward fashion. In particular, it enables us to easily draw 
plots of the relevant quantities for the various cases we 
will consider. 


III. MINIMUM FAILURE PROBABILITY FOR 
A FIXED DEGREE OF SEPARATION 

When the overlap of the final states is fixed, as in per¬ 
fect cloning, we argued above that a natural problem 
consists in deriving the minimum failure rate of the op¬ 
timal protocol, Qmin, as a function of one of the priors, 
say 771. In this section we address this problem by follow¬ 
ing the method employed in our derivation for cloning 
in m■ All the expressions below can be obtained from 
their analogs in [ 15 ] with the simple replacements s m —> s 
and s n —> s ', starting with the symmetric parametriza- 
tion of the curve Its lower half (for which q 2 < 91) 
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is parametrized as 


1 -xy- (-l)Vl - 3 ?a /1 - y 2 
2 


where 


1 — (1 + s')t 

x =- 77 -> y = 

s'/s 


* = 1 , 2 , ( 9 ) 

1 - (1 - s')t 


s'/s 


( 10 ) 


This parametrization arises from a change of variables 
that linearizes the unitarity constraint, which proved 
very convenient in m , where the advantages of its highly 
symmetric form were also apparent. The upper half of 
the curve @ can be obtained by applying the transfor¬ 
mation qi q2- However, without any loss of generality, 
we can assume that 0 < 771 < 1/2 (thus, 1/2 < 772 < 1 ), so 
only the lower half given by Eq. 0 can actually become 
tangent to the straight segment in Eq. 0 

Fig. [2] (a) shows plots of the unitarity curve [Eq. ^ 
plus the reflection 91 o 92] for s = 0.6 and s' = 0 . 05 , 
0 . 3 , 0.5 and 0 . 59 . For s' = 0 . 59 , very close to the 
value of s (small separation), the vertex of the curve ap¬ 
proaches the origin, which becomes a singular point in 
the limit s' —> s. As s' decreases (increasing separation), 
the curves approach the hyperbola 9192 = s 2 . It is ap¬ 
parent from the figure that the curves merge smoothly 
onto the lines q± = 1 and 92 = 1 for the larger values 
of s' . It becomes less obvious for small values of s' , such 
as s' = 0 . 05 . However a blowup of Fig. [2] (a) would reveal 
that this is so. A cusp at (s 2 ,l) and (l,s 2 ) arises only 
for s' — » 0 . 

It follows from our lemma, and it can be checked using 
Eq. that the slope of the lower half of the unitarity 
curve increases monotonically as we move away from the 
line qi = <72, where it has the value — 1 , and vanishes 
before we reach the line q± = 1 . The values of t at which 
the slope is —1 and 0 are, respectively, 


t-i 


1 — s'/s 1 — s' 2 /s 2 

1 — s' ’ 0 1 — s' 2 


( 11 ) 


So, there is a straight segment ([!]), with slope —*71/172, 
that is tangent to each point (91(f),(72(1)), t £ [t_i,t 0 ], 
of the unitarity curve parametrized by Eq. 0 - Since 
the slope of this curve is q' 2 (t)/q'^t ), where the prime 
stands for derivative with respect to t, a parametric ex¬ 
pression for 771 can be obtained from the equal slope 
condition —*71/172 = <7 2 (l)/<7i(l). The parametric expres¬ 
sion for Q m i n follows from imposing that (91 (1),<72(1)) 
must be a point of the straight segment ([I]), so Q m i„ = 
*7i9i(t) + *7292(1)- The final result can be cast as 


Vi 


/ 1 > Q min 

92 - Qi 


9291 - 9i92 
92 - 9 i 


l-i < 1 < lo, (12) 


where we have dropped the argument of qi(t) and q[{t) 
to simplify the equation. Further, one can check that the 
derivatives of qt(t) can be written as 


, = ^*(1 - qi) I 1 + 5' _ , 1 -s' \ 

i s'/s 1 


( 13 ) 


(a) (b) 




FIG. 2 . (a) Unitarity curves for different values of s'. 

The curves are symmetric under mirror reflexion along the 
(dotted) straight line 91 = 92, i.e., under the transfor¬ 
mation 91 92. (b) Minimum separation failure probabil¬ 

ity Q min vs. *71 (solid lines), for the same values of s' used 
in (a). In both figures the dashed lines correspond to full sep¬ 
aration/unambiguous discrimination ( s' = 0) and the value 
of the initial overlap is s = 0.6. 


Eq. ( fl 2 ] ) gives Q m in(* 7 i) in parametric form for 0 < s' < s. 
The solution for s' = 0 was already derived in the pre¬ 
vious section and for s' = s we have the trivial solu¬ 
tion Q min = 0 . These special cases can also be derived 


from Eq. ( 12 ) by carefully taking the corresponding lim¬ 
its. The values of Q m i n at the end points of this range 
follow by substituting t 0 and t_i, Eq. (Ill, into Eq. 0 . 
They are given by 


Q 0 — 92(^0) — 


1 - s' 


Q-i — 


1 -s'’ 


( 14 ) 


where Q m i n = Q-± holds for equal priors and Q min = Q 0 
for ?7i —► 0 (i.e., *72 -7 1). 

Fig. [2 (b) shows plots of the curves Q m in(* 7 i) for the 
same values of s and s' as the ones given above (solid 
lines). We see that Q m ; n is an increasing function of *71 in 
the given range [ 0 , 1 / 2 ], as one should expect. The figure 
also shows the failure rate for unambiguous discrimina¬ 
tion (dashed line), which coincides with Qmin for s' = 0. 
From the plots, it is clear that Q m i n is a decreasing func¬ 
tion of s ', again as it should be. 


IV. MAXIMUM SEPARATION 

In this section, we assume that *71, *72 are fixed given 
quantities and we focus on the relationships among the 
initial overlap, the final overlap and the maximum al¬ 
lowed failure rate. To find the explicit form of these 
relationships, we will need to develop a new geomet¬ 
ric view of both the unitarity constraint, Eq. © , and 
Q = *7i9i + *7292- We aim at a geometric representa¬ 
tion simple enough to grasp visually the solution and yet 
powerful enough to provide this solution analytically. We 
show below that the unitarity curve and the straight seg¬ 
ment of the previous sections can be mapped into conic 
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curves, in particular into families of parabolas and el¬ 
lipses respectively. This is arguably the simplest exten¬ 
sion to our geometric description of state separation. The 
desired transformation, similar in spirit to that in |[22j , is 
defined in terms of the new variables u and v as 


u = V?i92; 


v = 


gl + 92 

2 


(15) 


They are just the geometric and arithmetic means of the 
failure probabilities, q\ and q 2 - Under this transforma¬ 
tion the unitarity constraint becomes a parabola that can 
be conveniently written as 


1 + u 2 (u — s ) 2 
2 2s' 2 


(16) 


From this expression, one can immediately check that 
as s varies we obtain a family of parabolas whose envelope 
is yet another parabola, v = (1 + u 2 )/ 2, independently 
of s'. As s' decreases from its maximum value s' = s, 
the parabolas in Eq. (16) become thinner. For s' = 0 


they degenerate into the vertical segment u = s, 0 < v < 
(1 + s 2 )/2. These features are illustrated in Fig. | (a). 
Under the same transformation, Eq. (15), the line 


Q = mu + mm becomes an ellipse, which is most easily 
expressed parametrically in terms of the polar angle 9, 
measured relative to the axis v = 0 from the center of 
the ellipse. It is given by 


u = 


v = 


Q 

, cos 0, 

VI- A 2 

Q QA 


1 — A 2 1 — A 2 


sin (9, 


(17) 


where we have defined A = 772 — ^7i - It is clear from 
this expression that the eccentricity of the ellipse is only 
a function of the priors. For equal priors, A = 0, the 
ellipse degenerates into the horizontal segment v = Q, 
0 < u < Q, whereas for Q = 0 it collapses into the origin 
(u,v) = (0,0). As one increases Q, a family of similar 
ellipses is obtained. As they increase in size, their center 
moves up along the v axis. The line u = v is the envelope 
of this family, as one can easily check using Eq. (113- 
Fig. i (a) also illustrates these features. 

In terms of this conic geometry, optimality is again 
given by a tangency point, this time between ellipses 
and parabolas. Because of the features of these fam¬ 
ilies of conics, these points of tangency necessarily lie 
in the region between their envelopes, which is the gray 
area in Fig. i Fig. 1 (b) illustrates optimality. Given 
a maximum failure rate <3 max and some initial overlap s 
{Q max = 0.35 and s = 0.4 in the example considered in 
the figure), we plot the corresponding ellipse defined by 
Eq. © (dashed line). Among the various parabolas, 
characterized by the final overlap s' (the figure shows 
two of them, for s' = s and s' = s/2), the one that mini¬ 
mizes s' (solid line) has a unique point of tangency with 
the ellipse, thus giving us the solution, sU n . To keep the 
notation simple we will drop the subscript “min” wher¬ 
ever no confusion arises. 


(b) 




ry \ \s' — s 


\s’ = s/ 2 \ 


FIG. 3. (a) Unitarity parabolas, Eq. (16 1 , for different values 
of s, s' = s/7 (solid lines) and s' = s/1 4 (dotted lines). The 
dashed lines are the ellipses in Eq. ( |17[ ) for various values of 
the failure rate Q. The top boundary line to the gray region, 
given by v = (1 + u 2 )/ 2 , is the envelope of the solid and 
dotted parabolas. The bottom boundary line, i.e., the straight 
line v = u, is the envelope of the family of ellipses (dashed 
lines). The geometric solution to optimal separation falls in 
the gray region. In this figure 771 = 0.4. The degenerate 
curves for s' = 0 (dot-dashed vertical line) and A = 0 (dot- 
dashed horizontal line) are also shown, (b) Optimal (solid) 
and suboptimal (dotted) parabolas. The tangency point is 
also displayed. In this figure 771 = 0.3, s = 0.4 and Q = 
Qmax = 0.35. The optimal (minimum) value of s', which 
gives the solid parabola, turns out to be s' = 0.032. 


To find the condition that gives the tangency point, we 
first note that the slopes of the ellipse and the parabolas 
are given respectively by 


dv v A 

— = — =- . cot 0 , 

du u' yj\ — A 2 


dv u — s 

du U s' 2 


(18) 


where in the first line the primes stand for derivative 
with respect to the polar angle 9. The right hand side of 
these two equations must be equal at the tangency point. 
Moreover, the tangency point must belong to both the 
ellipse and the optimal parabola. Hence 


1+Asin# 1 Q 2 cos 2 9 1 / Qcos9 \ 2 

Q 1-A 2 = 2 + 2(l —A^^VyiUA^ 7’ 

A cot 9 1 — s' 2 Q cos 9 s 

y/1 - A 2 “ s' 2 y/i 3/P ~ V 2 ' 


where to obtain the first (second) equation we have sim¬ 
ply substituted Eq. ( fl7t into Eq. (16) [Eq. (18)]. Ide¬ 
ally, we would like to solve this system of equations by 
eliminating 9, which would lead to a closed expression 
relating s, s' and Q. Unfortunately, this involves solv¬ 
ing a high degree polynomial equation in cos 9. Instead, 
we look at it as a system of two equations with two un¬ 
knowns, s and s' (or Q and s') and keep 9 as a parameter 
describing the curve s'(s) [or .s'(Q)] in parametric form. 
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After some algebra, we obtain the simple expressions: 


can attain full separation (s' = 0). Past the critical value, 


y^l —Q ) 2 —(A+Qsinfl ) 2 
A + Qsind 


tan$, 


( 20 ) 


QA(l+sin 2 0) — (1 — A 2 — 2Q) sind 
S=- - -:-• ( 21 ) 


\J\ — A 2 (A+<2sin0) cosd 


The range of values of the parameter 9 in this equation 

in A ^ S' 


is — arcsin A < 9 < 0 max , where 


if Q < 1 - A, 


9 = ■ 

''max 


'max-’i . l-Q-A ( 22 ) 

arcsm-—- u Q > 1 — A. 


Q 


One can easily check the given minimum value of 9 by 
substituting in Eqs. (20 1 and (21) to obtain s' = s = 1, as 


it should be. Likewise, one can check that for 9 = d max 
one has s' = 0. The two cases in Eq. (22) reveal the 
appearance of the phase transition in the limit s' —> 0 
that we discussed in previous sections. If Q > 1 — A, 
substituting the second line of Eq. (22) in Eq. (19) we 
obtain s = [(2 Q + A — 1)/(1 + A)] 1 / . Solving for Q , we 
find that Q = 771 + s 2 r] 2 - This means that the condition 
Q > 1 — A is equivalent to 771 + s 2 r ]2 > 1 — A, which 
can be immediately seen to give 771 < s 2 /( 1 + s 2 ). So 
we obtain the second line in Eq. ([ 8 ]), corresponding to 
the “symmetry-broken phase”. If Q < 1 — A, namely, if 
s 2 /(l + s 2 ) < 771 , we have instead s = Q/\/l — A 2 . This 
equation can be written as Q = 2 ^ 771772 s. So, Eq. (22) 
has the same content as Eq. ([ 8 ]). Recall that we are 
assuming 771 < 1/2 < 1/(1+ s 2 ). The third line in Eq. ([s]) 
never applies under this assumption. 

Eqs. m an d ( pit are plotted in Fig. [4] (a) for two pos¬ 
sible priors: 771 = 0.1 (solid lines) and 771 = 0.5, i.e., for 
equal priors (dashed lines). From left to right, the max¬ 
imum allowed failure rate Q max is 0.2, 0.4, 0.6 and 0.8. 
We see that for small values of the initial overlap, s, one 


(a) 


(b) 




FIG. 4. (a) Plots of s' vs. s for 771 =0.1 (solid lines) and 
r/i = 0.5 (straight dashed lines) and for values of the failure 
rate. From left to right Q max = 0.2, 0.4, 0.6, 0.8. The dotted 
line is the (trivial) curve for Q ma . x = 0 , which is the straight 
line s' = s. (b) Minimum final overlap vs. maximum failure 
probability for various values of the initial overlap and the 
same two values of 771 used in (a). 


Scr 


Q 


max 





if Q < 2771 , 


if Q > 2771 , 


(23) 


full separation is no longer possible and s' increases (quite 
abruptly for small rfr). In the region s < s cr , the mar¬ 
gin < 5 max is not saturated, since the failure probability for 
unambiguous discrimination, Q ud, is smaller than Q max . 
For s > s cr we necessarily have to saturate the margin, 
i.e., Q = Q max- For equal priors (dashed lines) one can 
obtain the curves in explicit form from Eq. Q using that 
qi = <?2 = Q- 


s 


/ 


0 


s — 


Q 


max 


1 -Q 


max 


if S A f/max? 

if S > Qmax- 


(24) 


This expression could also be obtained by carefully tak¬ 
ing the limit A —> 0 in Eqs. (20) through (22). The figure 
clearly shows that separation becomes less demanding as 
we move away from the equal prior case. For Q max = 0, 
i.e., in the deterministic limit, we recover the trivial so¬ 
lution s' = s (dotted line). 


V. TRADEOFF BETWEEN MAXIMUM 
SEPARATION AND FAILURE RATE 


By solving the system Eq. (19) for Q and s' , we obtain 


a parametric expression for the tradeoff curve s'(Q ) in 
terms of the polar angle 9: 


Wl —A 2 


sin 9 


Q = 


^ A+sin 9 j 
\/l —A 2 (l+s 2 )cos 9— 2s(l + A sind) 
cos 9 

s\/l — A 2 + A s /2 cot 9 


(1 —s' 2 ) COS0 


(25) 

(26) 


Note that Eq. (25) is an expression for the square of the 


final overlap. To keep the formula for Q 1 Eq. (26), short, 
we use s' 2 as a shorthand for Eq. (251. The range of 9 in 
Eqs. (p5l) and (p6|) is: 


— arctan 


sA 


y/1- A 2 

where the upper limit of the interval can be written as 


0 


if Vi> 


^max * 


2s\/l - A 2 , r s 
-arccos^-—— , if 771 < 


l + s 2 ’ 
2 


(27) 


l-A + s 2 (l + A) 


1 + s 2 


The lower limit in the range of allowed 9 can be derived 
from Eqs. (25) and (26) by imposing that Q = 0 at s' = s. 













































Mirror 



FIG. 5. Six-port linear optics implementation of the optimal 
state separation protocol. The transmission (reflection) coef¬ 
ficients of the beamsplitters, BS1 and BS2 are given by the 
(off-)diagonal entries of the matrices in Eqs. (311 and (321, 
respectively. The input states are fed through ports 1 and 2 
as a superposition of zero and one photons in each port. The 
separated states are output through ports 1' and 2'. Port 3 
at the input is always in the vacuum state. A click in the 
photodetector placed in port 3' signals failure. 


The upper limit can be derived from Eq. (251 by impos¬ 
ing s' = 0. Once again, we see that a second order phase 
transition occurs in the limit of full separation: by sub¬ 
stituting the first (second) line of Eq. (27) in Eq. (26) we 
obtain Q = s\/\ — A 2 ( Q = [1 —A + s 2 (l + A)]/2), which 
is the first (second) case in Eq. ([ 8 ]). 

Fig. [4] (b) shows various plots of the separation 
VS. Qmax- As in Fig. i the plots are for ?y 1 = 0.1 
(solid lines) and for equal priors, 771 = 772 = 0.5 (dashed 
lines). For equal priors, there is the explicit formula for 
the curves given in Eq. (24). Again, we see that as 771 
gets smaller, departing from the equal prior value 1 / 2 , 
the states can be separated more for the same maximum 
rate of failure. As Q ma , x increases, the minimum over¬ 
lap gets smaller, as it should. When the margin Q max 
reaches the unambiguous discrimination value Qud we 


have s' = 0, attaining full separation. Larger values 
of Qmax are rather meaningless in this context, since they 
will never be saturated by an optimal protocol, which re¬ 
quires a failure rate of only Q = Q ud (< Qmax) to fully 
separate the input states. 


VI. A PHYSICAL IMPLEMENTATION: 

SINGLE-PHOTON MULTIPORT 
INTERFEROMETRY 

In this section we propose a physical implementa¬ 
tion of optimal state separation. The implementation is 
based on the dual-rail representation of qubits and single- 
photon multiport interferometry using only linear optics 
elements, namely, a mirror and two beam splitters, BS1 
and BS2. The measurements are carried out by three 
photodetectors. The setup is sketched in Fig.[5j 

The three input ports are labeled 1, 2, 3 in the fig¬ 
ure. A three-dimensional Hilbert space is spanned by 
the three orthogonal basis vectors corresponding to one 
photon in port i and vacuum in the other two ports. 
Thus, the basis vectors are |1) = a{|000) = |001), |2) = 
aJ,|000) = 1010) and |3) = ag|000) = 1100), where a] is the 
creation operator of the electromagnetic field in port i, 
i = 1,2,3 and |000) is the three-mode vacuum state. 
Similarly, for the output ports we have |1') = |001), 
12') = |010) and |3') = |100). 

In terms of these basis states, the input states are rep¬ 
resented as superpositions of |1) and |2). Note that the 
third port is always in the vacuum state at the input. 
Without loss of generality we choose the input states as 
|Vh) = |1) and (^2) = s|l) + a/ 1 — s 2 12), and the out¬ 
put states as |0() = |1') and \tp' 2 ) = s'11') + y/l — s' 2 12'). 
Then Eqs. ©> and © can be written as 


U |1) = VpI|1 , ) + v / 9T|3'), 

U (s|l) + \A-S 2 |2)) = yfc (s'11') + y/l - S' 2 |2'>) + V®|3'>, 


(28) 

(29) 


which corresponds to the choice |0)|a:o) = |3'). The de¬ 
tection of a photon in the output port 3' signals that 
separation failed. The state (02) can be produced in a 
standard way by sending a photon into a beam splitter 
with suitable transmission and reflection coefficients. 

For simplicity, we consider equal prior probabilities 
r/-| = r ]2 = 1/2, but the same setup can be used in the 
general case. As mentioned above, for equal priors we 
must have q\ = (72 = Q and p\ = P 2 = 1 — Q and the 


unitarity condition Eq. Q can be solved explicitly. The 
solution is given by Q = Q-i in Eq. ( fl4| . Substituting in 
Eqs. (28) and (29) we obtain two columns of the matrix 
of the unitary transformation U in the basis introduced 


above. The remaining column can be easily obtained im¬ 
posing unitarity. After some algebra we have 


[U] 



_ s—s' _ 

vdwhi+s) 


/ 


(l-M')(s-s') \ 

(W)(l+ S ) 


0 




(30) 


f as' ! (l-s)(s-s 7 y 

\V !- s ' V (1+«)(1-*') 


(l-s)(l-M') 

(l-s')(l+s) 




Using pill l!H] we can write U as the product U = M\ M-j. 
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where the matrices of M\ and M 2 are 




fv 

0 - 



[Mi] 

= 


0 1 

0 

> (31) 



w 

'fA » v 





/I 

0 



[M 2 ] 

= 

0 

/l+s' 

V l+s 

/ s—s' 

V l + s 

(32) 



1° 

s—s* 

V i+s v 

fm) 


We immediately 

recognize that the transformation Mi 


and M 2 can be implemented with beamsplitters, labeled 
in Fig. [5] by BS1 and BS2, respectively. The correspond¬ 
ing matrix elements provide the transmission (diagonal) 
and reflection (off-diagonal) coefficients of these beam¬ 
splitters. 

The degree of separation attained by the protocol can 
be certified by statistical analysis of the photon counts 
in the detectors placed in the ports V and 2', whereas 
those in the detector placed in port 3' provide the failure 
rate Q. 

Alternatively, one might consider the transformation 
provided by the set-up as a subroutine, probabilistically 
performing the requested state separation, as part of a 
larger protocol. One can achieve this by removing the 
detectors in V and 2' and feeding the output states into 
some subsequent unit for further processing. Hence, this 
implementation can be thought of as a separation module 
in a larger set-up. 


VII. CONCLUSIONS AND OUTLOOK 

In this paper we have addressed quantum state separa¬ 
tion for two known pure states with arbitrary prior prob¬ 
abilities. The degree of separation required by a prob¬ 
abilistic transformation determines its minimum failure 
rate. Thus, knowing the relationship between these quan¬ 
tities for arbitrary priors and arbitrary overlap of the in¬ 
put states is a valuable piece of knowledge for quantum 
information processing. It provides the ultimate limits on 
the processing of information allowed by nature and sets 
the performance scale for experimental implementations 
of such processing protocols. 

We have given a full account of state separation by fo¬ 
cusing separately on the various situations that one may 
encounter in quantum state processing. We first dealt 
with the optimization of protocols that have a fixed de¬ 


gree of separation, such as probabilistic perfect cloning. 
We have revisited, completed and extended our results 
in m- We have also given some technical details that 
were missing there. We have next considered the opti¬ 
mization of protocols for which a maximum allowed fail¬ 
ure rate, or margin, is given. We have computed the max¬ 
imum separation that a state transformation can possibly 
achieve as a function of the overlap of the input states 
and we have characterized the tradeoff between separa¬ 
tion and failure rate for fixed initial overlap. 

We have shown that a phenomenon analogous to a sec¬ 
ond order symmetry breaking phase transition arises in 
the limit of full separation, when the processed states be¬ 
come orthogonal. We have characterized it in the various 
situations discussed in the previous paragraph. Similar 
phase transitions have been discussed in connection with 
unambiguous discrimination of two or more states. The 
phenomenon arises from the high non-linearity of the uni- 
tarity constraints imposed by quantum mechanics. 

We have approached the optimization problems dis¬ 
cussed in this paper from a geometrical viewpoint that 
enabled us to gain a great deal of intuition about the so¬ 
lutions. This intuition has been the guiding line towards 
finding analytical results. Although a closed form for the 
solutions does not exist in the general case because of the 
high-degree non-linearity of the problem, our approach 
provides all the required relations between the relevant 
quantities in a clear and detailed way. The same geomet¬ 
rical approach has been applied in |15| and [21] where it 
proved equally powerful, and it can be applied to other 
optimization problems in quantum information process¬ 
ing where similar highly non-linear constraints arise. In 
this direction, we have some work in progress on prob¬ 
abilistic approximate cloning of two states and perfect 
cloning of three states. 
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